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Abstract 

The thermal conductivity of graphene nanoribbons (layer from 1 to 8 atomic planes) is investigated by using the nonequilibrium 
molecular dynamics method. We present that the room-temperature thermal conductivity decays monotonically with the number of the 
layers in few-layer graphene. The superiority of zigzag graphene in thermal conductivity is only available in high temperature region 
and disappears in multi-layer case. It is explained that the phonon spectral shrink in high frequency induces the change of thermal 
conductivity. It is also reported that single-layer graphene has better ballistic transport property than the multi-layer graphene. 



o 
o 



in 



X 



In the past decade, more and more attentions have been 
given to the question of what happens with thermal con- 
ductivity when goes to low-dimensional materials [1]. A 
two-dimensional materials-graphene (2D, in addition to its 
exceptional electric [3] and optical properties [HI], re- 
veals unique high thermal conductivity. Thermal conduc- 
tivity of single-layer graphene as well as of carbon nan- 
otubes is dependent on the chirality [6]. Recent theoretical 
studies suggest that the thermal conductivity of single -layer 
zigzag graphene is 20-50% larger than that of the single- 
layer armchair graphene [7]. However, whether the supe- 
rior thermal conductivity of zigzag graphene remains avail- 
able for multi-layer graphene has not got enough attention 
and concern. 

Additionally, experimental demonstrations have shown 
that the thermal conductivity gets a decrease at the two- 
to three-dimensional (2D to 3D) crossover of few-layer 
graphene [8J. The fact that the thermal conductivity of 
large enough graphene sheets should be higher than that of 
basal planes of bulk graphite was predicted theoretically by 
Klemens |9Q. Generally, thermal transport in conventional 
thin films still retains 'bulk' features because the cross- 
sections of these structures are measured in many atomic 
layers. Heat conduction in such nanostructures is domi- 
nated by extrinsic effects, for example, phonon-boundary 
or phonon-defect scattering [10]. A recent experimen- 
tal observation of high-quality few-layer graphene materi- 
als shows that the room-temperature thermal conductivity 
changes from ~2,800 to ~1,300 Wm _1 K _1 when the num- 
ber of atomic planes in few-layer graphene increases from 
2 to 4. It is explained that the observed evolution from two 
dimensions to bulk attributed to the cross-plane coupling of 
the low-energy phonons and changes in the phonon Umk- 
lapp scattering J8D. 

Recently, the method of molecular dynamics simula- 
tion has been successful in discovering thermal conduc- 
tivity and thermal rectification of the nanostructures 0]' 



I llll . This method, which builds the system from the bot- 
tom up, is useful to understand the intrinsic behavior, i.e., 
the phonon spectral behavior behind the significant change 
of a material's ability to conduct heat 1121 . In this pa- 
per, we will study the thermal conductivity of graphene 
ribbons (layer from 1 to 8 atomic planes) by using the 
nonequilibrium molecular dynamics method. By investi- 
gating the thermal conductivity and the phonon spectral 
behavior of multi-layer graphene, we try to reveal the cou- 
pling of the low-energy phonons from a more fundamental 
principle. The obtained results are of significance for un- 
derstanding the thermal properties in low-dimensional ma- 
terials and may open up few-layer graphene applications in 
the nanoscale thermal device such as thermal diode IU3I1 . 
thermal transistor IU4II and so on. 

In our simulations, we have used classical molecular dy- 
namics method based on the Tersoff-Brenner potential ltl5ll 
of C-C bonding interactions. The equations of motion for 
atoms in either the left or right Nose-Hoover thermostat are 
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where pi is the momentum and Ft is the force applied on 
the z-th atom. Q = 3A r /csT r 2 /2, where r is the re- 
laxation time, which is kept as lps. T is the dynamic 
parameter of the thermostat, T(t), which is defined as 
j^(vx(t) 2 + vy(t) 2 + vz(t) 2 ), where v(t) is the time- 
dependent velocity, is the instant temperature of the heat 
baths at time t. T (T L or T R ). The set temperature of the 
heat baths, are placed in the two ends of the graphene and 
the temperature difference is denoted as AT = T L — T R . 
For the convenience of comparison, here we set T L = 
350i^, T R = 300K, which are near the room tempera- 
ture. In order to avoid the spurious global rotation of the 
graphene in the simulation, we use fixed boundary condi- 
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tion in the two ends of the graphene. The fixed region and 
the heat baths occupy one layer and six layers of atoms re- 
spectively. N is the number of the atoms in the heat baths, 
ks is the Boltzmann constant and m is the mass of the car- 
bon atom. 

We integrate these equations of motion by a Verlet 
method 111 711 . The time step is 0.55 fs, and the simulation 
runs for 1 x 10 s time steps giving a total molecular dynam- 
ics time of 55 ns. Generally, T(t) can stabilize around the 
set value T after 2 ns for single-layer graphene and 20 ns 
for four-layer graphene. Time averaging of the tempera- 
ture and heat current is performed from 35 to 55 ns. The 
heat bath acts on the particle with a force —Ypi, thus the 
power of heat bath is —Tpf / m, which can also be regarded 
as the heat flux come out of the high temperature heat bath 
and injected into the low temperature heat bath. The total 
heat flux injected from the heat bath to the system can be 
obtained by J = J2 [- T Pi/ m i\ = -STNk B T(t), where 

i 

the subscript i runs over all the particles in the thermostat 
ifl^l . The final thermal conductivity is calculated from the 
well-known Fourier's law k = Jl/(ATwhn), in which I, 
w, and h (=0.335nm) are the length, width and thickness 
of each layer graphene, respectively, n is the number of the 
atomic planes. 

In Fig. 1, we explicitly observed the decrease of ther- 
mal conductivity as the number of layers changes from 1 
to 8. This change implies a crossover from 2D graphene 
to 3D graphite. For the single-layer graphene, the room- 
temperature thermal conductivity of zigzag nanoribbons is 
2276 Wm _1 K -1 , which is 43% larger than that of arm- 
chair nanoribbons. From the inset of Fig. 1, one can find 
that this attenuation relationship is not qualitatively varied 
with the graphene ribbon width. Another interesting phe- 
nomenon, that the zigzag graphene loses its superiority in 
thermal conductivity as the number of layers increases to 5, 
predicts that the chirality dependence of thermal conductiv- 
ity may disappear for 3D graphite. Here we have to point 
out that the superiority in thermal conductivity of zigzag 
nanoribbons is only available for the temperature larger 
than 220K. When the temperature is less than 220K, as dis- 
played in Fig. 2, we can observe that armchair graphene 
has better thermal conductivity than zigzag graphene. The 
temperature-dependent thermal conductivity of armchair 
and zigzag graphene can be interpreted by the number of 
phonons and the path for phonons passing. From the struc- 
ture of armchair and zigzag graphene, we can know that 
the path in armchair graphene (1.35/) is longer than that 
(1.15/) in zigzag graphene with the same length, where / 
is the length of graphene. Similarly, at the same width, 
the armchair graphene occupies more areas and particles 
than the zigzag graphene (the ratio is about 1.15). At high 
temperature, longer path means more phonons scattering. 
At low temperature, more particles mean more phonons to 
conduct heat. Therefore, the thermal conductivity of arm- 
chair is larger at low temperature and smaller at high tem- 



perature than that of the zigzag graphene. 
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FIG. 1: Layer-dependent thermal conductivity of few-layer 
graphene. The 1-layer nanoribbons refer to graphene and the 
5~8 layers nanoribbons are similar to utra-thin graphite. The in- 
set is the layer-dependent thermal conductivity for various width 
zigzag graphene. The average temperature is 325K. 
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FIG. 2: Temperature dependence of thermal conductivity of 
zigzag and armchair graphene for one and four layers atomic 
planes. 

Recent experiments suggest that thermal transport at the 
nanoscale is dominated by a ballistic rather than a diffu- 
sive mechanism lll8ll . For the single-layer graphene, fig- 
ure 3 illustrates that thermal conduction at low tempera- 
ture is dominated by the straight or bending mode, with a 
power law of ~T for armchair graphene and ~T 2 for zigzag 
graphene. The former is good agreement with the elastic- 
shell-based theoretical results in graphene with finite width 
ifl^l . The later shows that the chirality can change the 
power law relationship between the thermal conductivity 
and the temperature. We explained that the difference be- 
tween the power law thermal conductivity of armchair and 
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zigzag graphenes is mainly due to the different phonon 
scattering rates at the armchair and zigzag edges [7]. This 
power law relationship implies that the graphene conduct 
heat mainly through ballistic transport mode in this tem- 
perature region II 1911 . 

For the multi-layer graphene, the ballistic transport de- 
cays at the same temperature region and the thermal con- 
ductivity fluctuates near 500 Wm _1 K _1 . The crossover 
from 2D graphene to 3D graphite, which induces the de- 
crease of thermal conductivity of few-layer graphene, has 
been explained as the cross-plane coupling of the low- 
energy phonons and changes in the phonon Umklapp scat- 
tering [8]. For the single-layer graphene, i.e., in the ab- 
sence of cross-planes coupling, the thermal transport mode 
is of ballistic. Ballistic transport means less collision of 
phonons, however in the presence of cross-planes coupling, 
the phonons will scatters with the particles at the interface 
between the layers and then the thermal conductivity will 
decrease. Therefore, it is reasonable that the cross-plane 
coupling, which is similar to the inter-chain coupling in 
double Frenkel-Kontorova chains [20], will result in a sup- 
pression of thermal conductivity in graphene nanoribbons. 
Here, from the view of spatially-resolved phonon spectral 
behavior shown in Figs. 3(a) and (b), we present that the 
phonon occupies a wide frequency region for monolayer 
graphene in the absence of layer-to-layer coupling. How- 
ever, for four-layer graphene, in the presence of layer-to- 
layer coupling, as shown in Figs. 3(c) and (d) the phonons 
in high frequency region from 1,200 to 1,600 cm -1 shrinks 
and weakens at the edge of nanoribbons. In Refs. Il2ll 
and J2l|], the results about carbon nanotubes and graphenes 
have illustrated that the high frequency phonon modes, 
which can be affected by the edge sensitively, play a major 
role in the low-dimensional non-Fourier heat conduction. 
Obviously, the cross-planes coupling can induce unperfect 
and rough edges in the graphene. The edge roughness 
can suppress thermal conductivity strongly [12]. Thus, as 
the high frequency phonons shrinks with the layer-to-layer 
coupling, the material's ability to conduct heat will decay. 

In summary, the thermal conductivity of few-layer 
graphene strongly depends on the number of atomic planes. 
The reduction of the ability of few-layer graphene to con- 
duct heat is attributed to the crossover from 2D graphene 
to 3D graphite. We speculated this decreasing evolution in 
thermal conductivity is mainly due to the shrinking of high 
frequency phonon induced by the cross-layer coupling. We 
obtained the thermal conductance of zigzag and armchair 
nanoribbons as a power law relationship of the tempera- 
ture at low temperature. For the multi-layer graphene, the 
thermal conductivity is independent on the chirality. Al- 
though our calculations indeed show that the number of 
layer is disadvantageous to the thermal conductivity in few- 
layer graphene, the intriguing questions about the thermal 
conductivity of few-layer graphene are expected to stimu- 
late further experimental and theoretical investigations of 
phonon transport. 
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FIG. 3: Spatially-resolved phonon spectrum in the zigzag 
graphene. Edge (a) and middle (b) particle in the single-layer 
graphene, edge (c) and middle (d) particle in the four-layer 
graphene. The parameters are TV, — 350A, Tji — 300^. 
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